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HYPERNOVA AND GAMMA-RAY BURST REMNANTS AS TEV UNIDENTIFIED SOURCES 

KUNIHITO IOKA 1 , AND PETER MESZAROS 2 

ABSTRACT 

We investigate hypernova (hyper-energetic supernova) and gamma-ray burst (GRB) remnants in our Galaxy 
as TeV gamma-ray sources, particularly in the role of potential TeV unidentified sources, which have no clear 
counterpart at other wavelengths. We show that the observed bright sources in the TeV sky could be dominated 
by GRB/hypernova remnants, even though they are fewer than supernova remnants (SNRs). If this is the case, 
TeV SNRs are more extended (and more numerous) than deduced from current observations. In keeping with 
their role as cosmic ray accelerators, we discuss hadronic gamma-ray emission from ir° decay, from (3 decay 
followed by inverse Compton emission, and propose a third, novel process of TeV gamma-ray emission arising 
from the decay of accelerated radioactive isotopes such as 56 Co entrained by relativistic or semi-relativistic jets 
in GRBs/hypernovae. We discuss the relevant observational signatures which could discriminate between these 
three mechanisms. 

Subject headings: cosmic rays — gamma rays: bursts — gamma rays: theory — radiation mechanism: non- 
thermal — supernova remnants 



1. INTRODUCTION 

Observations at a new wavelengths have always led in as- 
tronomy to new discoveries, such as radio pulsars, gamma- 
ray bursts (GRBs), etc. The TeV gamma-ray sky is likely to 
yield similar surprise s, since among the rapid ly growing num- 
ber of TeV sources dAharonian et al.l 12006) the most abun- 
dant category is of a mysterious nature, the so-called TeV 
unidentified sources (TeV unl Ds), which have no clea r coun- 
terpart at other wavelengths ( Aharoni an et alj 12005. 2008; 
Mukherjee and Halpern 2005). Our knowledge of TeV unlDs 
is very limited: 

(1) So farA/ un iD ~ 10-30 TeV unlDs have been observed. 

(2) They generally lie close to the Galactic plane, suggest- 
ing a Galactic origin. 

(3) They are extended, Ail ~ 0.05-0.3°. 
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implying a rarer type of source. 

Long GRBs are actually rare SNe endowed with rel- 
ativistic je ts, which also ar e expected to leave SNR-like 
remnants dPerna et al.l l2000t lAval & Piranl 12001b . Even 
so, the expected number of GRB remnants is so small 
that previous works have addressed only specific sources 
(Iloka etal.ll2004t lAtoyan et all 120061) . However, recent ob- 
servations suggest that hyper-energetic SNe, the so-called 
hypernovae which are sometimes associated with GRBs 
such as SN1998bw/GRB980 425, SN2003dh/G RB030329 
and SN20031w/GRB031203 (Maedaand NomotoS l200l). oc - 
cur more frequently than GRBs (Guetta & Delia Valle 2007). 
The hypernova rate may be even higher if we are miss- 
ing GRB-unassociated hypernovae like SN1997ef. On the 
other hand, a larger fraction of SNe may be endowed with 
slower or semi-relativistic jets, causing only low-luminosity 



(5) They have a power-law spectrum with index of 2.1-2.5. (LL) GRBs. Actually, a significant energy — 10 erg 



(6) Some TeV unlDs have strong upper limits in X-rays 
with a T eV to X-ray flux r atio o f Fj e y/Fx > 50 from 
Suzaku dMatsumoto et all 120071: iBamba et all [ 2007) 
and in radio with Fj e y / F ms n ~ 10 3 dAtovan et alj 2006; 
ITian et al.ll2008h . 

The high energy nature of TeV unlDs naturally leads us to 
consider the cosmic ray (CR) accelerators as their possible 
origin. Since the dominant CR sources are supernova rem- 
nants (SNRs), TeV unlDs may be related to SNRs, in partic- 
ular old ones, which are expected to be less luminous at other 
wavelengths ( Yama zaki et alJl2006l) . However, the required 
galactic energy budget of TeV unlDs is only 
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was released in the form of a mildly relativistic ejecta in 
the very faint GRB980425/SN1998bw. The recently dis- 
covered LL GRB060218/SN2006aj coul d be also be pro- 
duced by a slower o r semi-relativistic jet dToma et al.ll2007t 
IWaxman et al.l 120071). occurring at a — 10 times higher rate 
than GRBs dGuetta & Delia Valid 120071: iLiang et all 120071: 
Soderber g et al.H2006d) . Such semi-relativistic jets could be 
intermediate between choked and breaking-out jets emerging 
through the progenitor's outer envelope. 

An interesting point in TeV unlDs is that their emission 
may be related with hadronic processes, since a simple lep- 
tonic process for TeV gamma-rays by inverse Compton (IC) 
would require too much synchrotron emission than what is 
observed in TeV unlDs. Therefore we co nside r the hadronic 
processes associated with the ir° decay (§ 12. U f or th e hyper- 
nova shocks, and with the (3 decay mechanism (§ 13. 11 1, as well 
as a ne w mechanism of accelerated radio-isotope (RI) decay 
(§ I3.2| | for the GRB jets. In the RI decay model, TeV gamma- 
rays are produced by the Lorentz-boosting of the MeV decay 
gamma-rays of accelerated RI, such as 56 Co entrained by the 
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jets. To our knowledge, this is the first discussion of the decay 
of accelerated RI as a mechanism for generating Te V gamma- 
rays in astrophysical sources. In this case, the source of the 
gamma-ray energy is injected initially at the base of the jet. 
There is no need for target matter or photons, as in the usual 
processes such a s the ir° decay, inverse C ompton (IC) or pho- 
todisintegration (lAnchordoqui et aDl2007h . 

In this paper, we investigate the high energy implications of 
GRB/hypernova remnants for TeV sources, in particular for 
TeV unlDs. Even though the total number of GRB/hypernova 
remnants is less than that of SN Rs, the observed number in 
the TeV sky could be larger (§ 12. U . This in turn predicts 
more extended (and more numerous) TeV SNRs than cur- 
rently observed, which may be discovered by expanding the 
search region even with the current instruments. We use the 
units kg = h = 1 and Q x = 2/10* in cgs units unless otherwise 
stated. 

2. HYPERNOVA SHOCKS 
2.1. 7T° decay model 

The simplest and most plausible process for TeV gamma- 
ray emission associated with CRs is the 7r° decay from pp 
interactions between the interstellar medium (ISM) and CRs 
accelerated by th e SN/hypernova shocks . Although the 7r° 
decay for SNRs (Naito & Takahar a1ll994t iDrurv et al.lll994l) 
and GRB remnants ( Ato van et al.l 120061) has been discussed 
in detail, there are few works dealing with tt° decay in hyper- 
nova remnants J _excepjj[ or the specific source HESS J1303- 
631 (lAtoyan etal.l l2006). Since hypernovae are more ener- 
getic than SNe and GRBs and more frequent than GRBs, it is 
worth considering their implications for TeV unlDs. 

By scaling the SNR calculation dNaito & Takaharal fl994t 
IDrurv et al.lll994h . we obtain a TeV gamma-ray flux 

£l F Ei ~ 10- 12 C-i£ 5 i«</rok P c erg s" 1 cm" 2 , (3) 

comparable to the observed values. Here n is the ISM density 
in cm" 3 , ofiokpc = d/lOkpc and £ is the fraction of the total CR 
energy E per logarithmic energy interval. Since the kinetic 
energy of hypernovae is huge, ~ 10 52 erg, a CR energy of 
E ~ 10 51 erg is reasonable, for a conventional CR acceleration 
efficiency of ~ 10%. 

Since the flux in equation (0 is independent of time, most 
sources are likely old, e.g., f age ~ 10 5 yr old. Such old SNRs 
are possible TeV unlDs because the maximum energy of pri- 
mary electrons is so small that the leptonic emission is sup- 
pressed. As shown in Fig. [T] the synchrotron emission from 
secondary electrons can also be belo w upper limits for TeV 
unlDs in the radi o to X-ray bands (lYamazaki et al.1 120061: 
lAtoyan et al.ll2006h . 

Remarkably, the observed number of old hypernova rem- 
nants may be comparable to that of the more numerous old 
SNRs, if the observations are flux-limited. For old remnants, 
which are larger than the angular resolution [see equation (0], 
the source size should be taken into account. As the search re- 
gion is expanded, more background is included. The sensitiv- 
ity is proportional to the inverse square of the background, for 
background-dominated counting statistics. Therefore, the flux 
sensitivity to an extended source with a physical radius r is 
given by f e extend = F£° mt (r/d9 cut ), where Fi° mt is the sensitiv- 
ity to a point source and 9 cut is the angular c ut in the analysis 
dKonopelko et al.ll2002t iLessard et al.H200H) . With equation 



© and 

F '° m idt) <F ^ End ~ 2 > (4) 

we have the maximum distance to a source as 

rfmax oc E n r~ l (5) 
and hence the observable volume in the Galactic disk is 

V oc c/ 2 ax oc E 2 n 2 r- 2 oc £ 8/ V 2/s , (6) 

where we use r oc E^rT 1 ^^, which is approximately cor- 
rect even for the radiative phase of the remnants. Since the 
hypernova energy is ~ 10 times larger, the observable num- 
ber of hypernova remnants is larger than that of SNRs by 

^hnrVhnr lO-V 1 ■ 10 8 / 5 Q4 

where R SNR ~ 10" 2 yr _1 and R HNR ~ lO^yr -1 (~ 7% of the 
SNe Ibc ra te) are the event rates of SN e and hypernovae, re- 
spectively dGuetta & Delia Vallel2007l) . and we assume a sim- 
ilar age f age and ISM density n for both remnants. The actual 
numbers in our Galaxy would be /^sNR^age ~ 10 3 (SNRs) and 
^HNR^age ~ 10 (hypernova remnants) for f age ~ 10 5 yr. Since 
the observed number of TeV unlDs is also ~ 10, we may be 
reaching the farthest hypernova remnants in our Galaxy. 
The angular size of the f age ~ 10 5 yr old remnants, 

L„mz„w(-±X t (8 ) 

d d V 10k P c / 

is consistent with TeV unlDs for hypernovae, 3 but the same 
may not be true for SNe since the observable distance d max is 
smaller and hence the angular size is more extended [r/d ~ 
2°(lkpc/ii)]. Such an extended object can be found only if 
the search region is expanded to the degree scale in the analy- 
ses. Therefore, if TeV unlDs are hypernova remnants, we can 
predict more extended (and more numerous) TeV SNRs than 
observed, which may be discovered by expanding the search 
region even with the current instruments. 

The flux of hypernova remnants could be further enhanced 
if the density n is higher around hypernovae than around 
usual SNe. The GRBs and SNe Ic, which are the only 
type of core collapse SNe associated with GRBs and hyper- 
novae, are far more concentrated in the very brig htest re- 
gions of their host galaxies than are the ohter SNe ( Kellv et al. 
120071; iFruchter et al]|2006l) . suggesting that GRBs and hyper- 
novae are associated with star forming regions with high den- 
sity. Conversely, hypernovae could be runaway massive stars 
ejected several hundreds of parsec away from high density 
stellar clusters (Hammer et al. 2006). In this case, no flux en- 
hancement is expected. 

It is necessary to check any possible side-effects of this 
mechanism at other wavelengths which might violate the 
hypthesis that they remain unidentified. In particular we need 
to check what is the possible X-ray emission from electrons 
in these sources. We may constrain the electron accelera- 
tion in the old hypernova remnants. In Fig. [T] we have plot- 
ted the bremsstrahlung, synchrotron, and IC emission from 
primary electrons with the same spectral index p = 2.2, for 

The CR diffusion could increase the source size JAtovan et aT]|200a) . 
although the diffusion coefficient has large uncertainties. 
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FIG. 1. — Flux from the n° decay via pp interactions between the ISM 
and CRs accelerated by the hypernova shocks, compared with the Fermi and 
HESS sensitivities. We assume a remnant of age ( a gc = 10 5 yr at d = 10 
kpc with a CR energy E = 3 X 10 erg and the CR spectral index p = 2.2 
in the energy range m p c 2 < e p < WprripC 1 . We also show the synchrotron 
emission from n° decay positrons and electrons for B = 3/^G, compared with 
the observational upper limits for an X-ray to TeV flux ratio of 10~ 2 and a 
radio to TeV flux ratio of 10~ 3 . We use a code of Kamae et al. 12006) for 
calculating pp interactions. In addition, we also plot the bremsstrahlung, 
synchrotron, and IC emission from primary electrons with the spectral index 
p = 2.2 and the electron-to-proton ratio of (ej p)ioGeV ~ m e /m p at lOGeV. 

an electron-to-proton ratio of (e//?)ioGeV ~ m e /m p ~ 10~ 3 at 
lOGeV. The radio observations limit the (e/p)ioGeV ratio to 
< 10~ 3 , which is somewhat smaller than the observations 
~ 0.02. Note however that the injection of thermal electrons 
into the acceleration process is poorly understood, and the 
observed electron CRs might be produced by other sources 
such as usual SNRs or pulsars. In any case, the leptonic 
model for TeV unlDs seems unlikely. Note that we take the 
Coulo mb losses for low energ y electrons into account accord- 
ing to lUchivama et all d2002l) (see also iBaring et all (119991) ) 
and this makes the bremsstrahlung spectrum hard in the X- 
ray band. On the other hand, the thermal bremsstrahlung 
emission is not bright in the X-ray band, in contrast with 
the usual young SNRs dKatz & Wa xman 2007), because the 
temperature of old SNRs is much below the X-ray energy, 

T ~ 0.04 keV ^ 2/ 5 5 2 n- 2 / 5 (f/10 5 yr)" 6/5 (lReynoldsll2008l) . 

3. GRBJETS 

3.1. [3 decay model 

A fraction of SNe/hypernovae is associated with long GRB 
and their jets. Another process for generating TeV gamma- 
rays is the (time-delayed) (3 decay of the neutron compo- 
nent of the CR outflow accel erated by the j ets followed by 
inverse Compton scattering (lloka et al.l 120041) . We have ap- 
plied this process to the S NR W49B (G43.3-0 . 2), which is 
a possible G RB remnant dKeohane et al.1 120061 iMiceli et al.1 
|200g[2008). Since the (3 decay can occur outside the rem- 
nant, the jet-like emission would appear outside of the SNR. 
Note that the GRB-associated SNe are not necessarily ener- 
getic hypernovae since the luminosity distribution of GRB- 
associate d SNe is statistically c onsistent with that of local 
SNe Ibc (ISoderberg et al.ll2006bl) . 

As shown in Fig. [2] old jet remnants with f age ~ 10 5 yr are 
potentially TeV unlDs since the usual SNR emission goes 
down at this age, while electrons that emit ^TeV gamma- 
rays (i.e., 7 e ~ 10 7 el ectrons) need f C Q i ~ 10 s yr to cool 
(see Eq. (3) and (4) in lloka et al .1120041) . (However a spec- 
tral cutoff at ~ 10 TeV may be seen at this age.) Then the 



FIG. 2. — Flux from the IC scattering of CMB photons by the /3 decay 
electrons via the CR neutron component in the GRB jets, compared with the 
Fermi and HESS sensitivities. We assume a remnant of age ( a ge = 10 5 yr at 
d = 10 kpc with a CR energy E = 3 X 10 51 erg [i.e., an old remnant version 
of Model (I) in lloka et all tjoOli l and a CR Lorentz factor 10 s < 7,, < 10 9 . 
We also show the synchrotron emission from (i decay electrons for B = 3/jG, 
compared with observational upper limits for an X-ray to TeV flux ratio of 
10~ 2 and a radio to TeV flux ratio of 10~ 3 . 

expected total number is 0.1-1 for a Galactic GRB rate of 
- lO^-lO" 6 yr" 1 dGuetta & Delia Valid l2Q07h . The actual 
number would be larger because a larger fraction of SNe may 
have lower-lumin osity (LL) jets that were not identified (e.g., 
Toma et al. 2007), and LL jets can also accelerate CRs to en- 
ergy ~ 10 16 " 17 eV (e.g.. lMurase"eTaIll2006l [2008h . which is 
necessary for the TeV gamma-ray emission. The most op- 
timistic rate consistent with the late-time ra dio observations 
is ~ 1 0" 4 yr -1 (- 10% of the SN Ibc rate) (ISoderberg et alj 
2006a), yielding a total number of ~ 10, comparable to that 
of TeV unlDs. All of these remnants could be detected, i.e., 
<^max ~ lOkpc, if the total energy in accelerated CRs is com- 
parable to the typical GRB energy. 

3.2. Radio-isotope (RI) decay model 

The chemical composition of jets which may be associated 
with SNe/hypernovae is currently unknown. One interest- 
ing possibility is that they entrain radioactive isotopes (RI), 
in particular 56 Ni and 56 Co. The SN/hypernova shock leads 
to explosive nucleosynthe sis, predominantly forming 56 Ni via 
complete silicon burning ( Ma eda and No moto 2003). The in- 
ner 56 Ni would fall back onto the central engine (in those ob- 
jects where there is one), and a fraction of accreted 56 Ni could 
be ejected with the ensuing jet. The photodisintegration of 
heavy nuclei would be suppressed in the LL jets whose tem- 
perature is less than ^MeV. Even if heavy nuclei may disinte- 
grate in the innermost accretion disk, 56 Ni could be produced 
again in a cooling wind (MacFadven 2003). 56 Ni could be 
also entrained into the jet from the surroundings. The subse- 
quent photodisin tegration can be avoided in large fract ion of 
parameter space dWang et al.l l2008: Mur ase et al.ll20()8). The 
spectrum and light curve modeling of GRB-associated hyper- 
novae suggests that the total amount of 56 Ni is larger than that 
in ordinary SNe and h eavy elements are aspheric ally ejected 
along the jet direction (Maeda and Nomoto 2002). 

56 Ni decays int o 56 Co in t^\ln2 = 6.1 days essentially by 
electron capture dMcLaughlin & Wijersl 120031) . The decay 
proceeds primarily through the excited states of 56 Co, emit- 
ting ~ 2 MeV gamma-rays. 56 Co also decays primarily by 
electron capture into excited 56 Fe with a half-life of fc ln2 = 
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77.2 days, 

56 Co + e--> 56 Fe* + i/ e , (9) 
followed by e RI ~ 2 MeV gamma-ray emission, 

56 Fe *_ > 56p e + 7 (1Q) 

If the nuclei are fully ionized, they decay by emission of a 
positron (j3 + decay), e.g., 

56 Co^ 56 Fe* + e + +^, (11) 

followed by gamma-ray emission, since free electron capture 
is negligible in our case. The half-life for ionized 56 Co is a 
factor of ~ 5 higher than for 56 Co with electrons, while that 
of ionized 56 Ni is longer than 3 x 10 4 yr. Since nuclei ca n 
recombine in an expanding jet (McLaughlin & Wiiers 2003), 
the decay chain 56 Ni — > 56 Co — > 56 Fe can start from the 56 Ni 
electron capture. 

The outflowing jet is shocked as it interacts with the ISM, 
leading to particle acceleration. Unlike in the ordinary SNe, 
the RI acceleration can occur before the RI decay in the 
objects considered here, involving jets reaching relativistic 
speeds. Our interest is in the 56 Co acceleration, since the ac- 
celerated nuclei are completely ionized and the ionized 56 Ni 
is essentially stable (the half-life is longer than 3 x 10 4 yr). 
For the 56 Co acceleration, the shock radius should be smaller 
than the 56 Co decay length, r < cfc (ln2)/3r and larger than 
the 56 Ni one, r > ctmQn2)f3T, where T is the Lorentz fac- 
tor of the jet. Since the shock radius may be estimated by 
Ej ~ (4Tr/3)r 3 nm p c 2 p 2 T 2 , we have 

2E 1 J 5 \n- l/5 <f3T< 8£ / 1/ 5 5 1 «- I/5 (12) 

where Ej is the total energy of the jet. The Lorentz factor 
ranges between that of GRBs and hypernovae. Note that the 
shock radius is re latively large > 10 16 cm that no heavy nuclei 
may disintegrate dMurase et aljiooalXVang et al.ll2008l) . 

Once 56 Co is accelerated, the observed mean lifetime is ex- 
tended by the Lorentz factor, 

h = 7Aco ~ 10 5 yr75 (13) 

enabling long-lasting emission, where f;c = 5fc is the half- 
life of ionized 56 Co. The observed energy of the decay 
gamma-rays is also Lorentz-boosted to 

£ 7 = 7£Ri~0.2TeV 75 , (14) 

potentially applicable to TeV unlDs. The ratio of gamma- 
ray energy to the 56 Co CR energy is about / = Sri / 56m p c 2 ~ 
4 x 10~ 5 . The gamma-ray flux is then estimated as 

= ~ 10- I2 erg s" 1 cm" 2 ^, (15) 

interestingly comparable to that of TeV unlDs, where a^pc = 
c//3kpc and E is the total energy of 56 Co CRs. In a sense, this 
is a relativistic version of the SN emission since the energy 
source is the RI decay. 

In the RI model, a similar way to the j3 decay model, 
the total source num ber is 0.1-1 for a Galactic GRB rate 
of - 10" 5 -10- 6 yr" 1 dGuetta & Delia Valid l2007h . while the 
most optimistic number is ~ 10 for a rate consistent with 
the late-time radio observations ~ 10~ 4 yr -1 (~ 10% of the 



SN Ibc rate) dSoderberg et al.ll2006ah . We can detect a frac- 
tion ~ (3kpc/10kpc) 2 ~ 0.1 of these sources, if CRs comprise 
an energy comparable to the typical GRB energy. For a rare 
event with high energy and/or with reacceleration by the hy- 
pernova shock, the observable distance may be farther. 

We note that other RI species may also contribute. In par- 
ticular, the 57 Ni emission via (3 + decay dominates in the young 
(~ 10 3 yr) remnants because its half-life 35.6 hr is ~ 10 2 
times shorter than 56 Ni though its yield is less by 0.1-0.01 
dMaeda and Nomotoll2002l 120031) . 

The shocked 56 Co is accelerated to a power-law spectrum 
dNco oc e^deco for eco < Eco.max- The maximum energy can 
be e C o,max ~ 3 x 10 17 eV (e.g jMurase et al.l2008l) . which pro- 
vides e 7 , max ~ / x (3 x 10 17 eV) ~ 10 TeV decay gamma-ray. 
Hereafter we adopt this maximum energy. The Larmor ra- 
dius of 56 Co CRs is ~ 4(e C o/10 17 eV)B:^ pc, leading to the 
isotropic emission of decay gamma-rays. 

As shown in Fig. [3] the decay gamma-rays have a unique 
spectrum, 

£l F Ei oc e- p+1 exp (-^j (16) 
with an exponential cutoff 4 below 

e - £ «(i)- 02TeV fe) (17) 

and a power-law above it, where the power-law index is softer 
than that of parent CRs by one. This is because the low en- 
ergy CRs have already decayed while the high energy CRs are 
decaying with a rate (d/dt)exp(-t /jhc ) oc e~ l exp(-e p /e 7 ). 
The spectrum of TeV unlDs is characterize d as a power- 
law with index 2.1-2.5 dAharonian et al.ll2008l) . which corre- 
sponds to p = 1.1-1.5 in equation ( fT6l ). This is slightly harder 
than the p = 2 expected in the test particle acceleration case. 
The harder spectrum could suggest non-linea r shock effect s 
in jets where the CR pressure is nonnegligible (Ellison 2001). 
This is consistent with the fact that the necessary CR energy 
is comparable to the total energy of jets, E ~ Ej. 

In Figure [3] we also show the synchrotron and IC emis- 
sion from decay positrons in equation ( fTTI ), which have an en- 
ergy comparable to that of decay gamma-rays. The positron 
spectrum is dN e oc e~ p de e for Tm e c 2 < e e < e p and dN e oc 
e'^dEe for e p < e e < e max - The synchrotron frequency is 
= qBe 2 e /2Tnnlc 5 ~ lO^B-efe/O.lTeV) 2 eV while the IC 
frequency is v ic ~ 9e CMB {e e /m e c 2 ) 2 - 10 8 (e e /0.1TeV) 2 eV 
where we consider c osmic microwa ve background (CMB) as 
main target photons dloka et alj|2004l) . 

However the cooling time t c = 3m 2 c 3 /4arS e U ~ 
10 7 (e e ,/0.1TeV)~ 1 t/Ti 1 2 yr is usually longer than the 
decay time f 7 in equation dTJb for e e < e max ~ 10 
TeV where U is the total energy density of magnetic 
fields and CMB. Then the luminosity is suppressed by 
tj/t c — 10" 2 (f 7 /10 5 yr)(e e AXlTeV)£/_i2 compared with decay 
gamma-rays in equation ( U5l ). This is favorable for the inter- 
pretation of such sources as TeV unlDs. Note that the ^eV 
photons by decay positrons come from an extended region in 
which the optical background dominates the remnant flux. 

4 The actual cutoff is a power low s~/F £y oc e 7 for e 7 < e p because 
nuclei also make off-ax is emission outside the relativistic beaming angle 
I Ioka & Nakamura 2001). This does not change our results. 
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FIG. 3. — - Fl u x of the deca y gamma-rays from accelerated Rl 56 Co CRs in 
equations UOt . )15t and 1161 . compared with the Fermi and HESS sensitivi- 
ties. We assume a remnant of age / age = 10 5 yr at d = 3 kpc with a CR energy 
E = 3 X 10 51 erg and spectral index p = 1.5. We also sho w th e synchrotron 
and IC emission from decay positrons given by equation (TTJ for B = 3/xG 
and the cosmic microwave background, compared with observational upper 
limits for an X-ray to TeV flux ratio of 10~~~ and a radio to TeV flux ratio of 

irr 3 . 

The accelerated Rl CRs may suffer the adiabatic 
energy losses during th e expansion of the remnant 
(Rachen & Meszarosl ll998l) . which will greatly suppress the 
Rl emission. This is a potential problem for the Rl decay 
model. However the CR escape from a remnant is not fully 
understood yet. For instance, it is common to assume the 
free escape of CRs in models where the GRB reverse shocks 
are the sites for the ultra hi gh en ergy CR production (e.g., 
IWaxmanl20 06; Murase_etalJ[2008). In the same sense, the Rl 
decay model is a viable possibility for TeV unlDs. 

4. DISCUSSION 

We have discussed the TeV gamma-ray emission from the 
7T° decay, (3 decay and the radio-isotope (Rl) decay mech- 
anisms in GRB and/or hypernova remnants, in their possi- 
ble role as TeV unlDs. There is evidence that typical GRBs 
predo minantly occur in galaxies with less m etals than our 
own dFruchteretal.1 120061: IStanek et afl 120061) . although the 
correlation between t he host metallicity an d the GRB en- 
ergy is not confirmed (Savagl io et al.ll2008allbl) . On the other 
hand broad-lined SN e Ic also inhabit more metal-rich galaxies 
(Modiaz et al. 2008), suggesting that at least hypernovae and 
possibly GRB/SNe with LL jets can occur in our Galaxy. In 
addition, the evidence that GRBs occur in low metal regions 
could be biased because the GRBs in the metal rich region 
tend to have no optical afterglows due to the dust absorption 
(i.e., dark GRBs, which is about half of all GRBs). Then it be- 
comes difficult to identify the hosts for measuring the metal 
abundance. 

To discriminate amongst these models, the imaging of the 
gamma-ray morphology would be useful. The 7r° decay 
model predicts a shell structure, since the density inside the 
remnants is low, while the (3 decay model predicts an elon- 
gated structure, and the Rl decay model predicts a center- 
filled structure. For remnant ages f age < 10 5 yr, the SNR 
may be observed also at other wavelengths. HESS J1731-347 
may be such an example, which is likely to be an old SNR 
with a large gam ma-ray to radio flux ratio Fr e v /-Fradio ~ 33 
dTian et al.ll2008b . HESS J 1834-087 may be also an old rem- 
nant with an elongated gamma-ray emission outside the SNR 
dAharonian et alj 120061) . which could be i nterpreted in the 
conetxt of a /3-decay model (but see also IMukherie e et al.l 



(2008)). 

The 7T° decay model predicts detectable GeV emission if 
the spectral index of 2.1-2.4 continues down to the GeV re- 
gion (see Fig. 0. Other m odels predict relatively low lev- 
els of GeV emission (Fig. [3j Ioka et al. 2004; Yamazaki et al. 
l2006tlAtoyan et al.l20 06). although rare young remnants may 
be bright if e p ~ 20GeV(f/10 4 yr) in the Rl model. This may 
be tested by the Fermi satellite. Note that the GRB prompt and 
afterglow emission will mask the very early Rl decay emis- 
sion, which becomes prominent later. 

The implied TeV neutrino flux (~ gamma-ray flux) is not 
enough for detection by current facilities. However, in the fu- 
ture one may in principle test models through the flavor ratio. 
We expect ir + — ► fi + + v ^ — > e + + v e + + in the 7r decay 
model, no neutrinos in the (3 decay model (since neutrinos are 
beamed and usually off-axis), and 56 Co — > 56 Fe* + e + + v e in 
the Rl decay model. 

As discussed at the end of § 12.11 the radio observations limit 
the (e//?)inGeV ratio to < 10~ 3 in the ir° decay model. For 
the Rl decay model, the CRs need to freely escape from the 
acceleration site befo re suffering the adiabatic energy losses 
(see the end of § 13.21 ). 

The implied CR energy budget is less than ~ 10% ~ 
(10 times energy) x (10~ 2 times rate) of the standard SN CR 
energy budget in our models. The CRs above the knee at 
~ 3 x 10 15 eV, however, could be produce d mainly by extra- 
galactic and/or Galactic GRBs/hyperno vae ( Wick et alll2004t 
iMurase et al.l 120061: IWang et al.l 120071: iBudnik et al.l 120081) . 
The chemical composition is increasingly richer in heavy 
nuclei above th e knee ~ 3 x 10 15 eV to the second knee 
- 6 x 10 17 eV (lAntoni et all 120051: lAbbasi et al J 120051) . and 
possibly above ~3 x 10 19 eV (Unger et al. 2007). These may 
be accelerated by jets as in the Rl decay model. 

The GRB remnants may be also responsible for the ex- 
cesses of cosmic -ray positrons and electrons recently ob- 
served by t he PAMELA and ATIC/PPB-BETS experiments 
(Ilokal2 008: Fuji ta et alJ2009l) . Since the electron and positron 
sources should be nearby (less than ~ lkpc away), it may be 
difficult to detect these sources by gamma-rays. 

Recently, wide-field optical surveys have been discovering 
exceptionally luminous SNe, such as SN 2 005ap, SN 20 08am, 
SN 2006gy, SN 2006tf and SN 2008es dMiller et al.ll25c>8h . 
The SN 2008es yields a radiated energy of > 10 51 erg and 
possibly a total energy of ~ 10 52 erg, comparable to the hyper- 
nova energy. The rate could be also comparable to the hyper- 
novae rate although only an upper limit to the rate is obtained 
(Miller et al. 2008). Therefore such luminous SNe could also 
become TeV unlDs like the hypernova case, although the lu- 
minous SNe are type II and may not be related with GRBs. 

In conclusion, we propose that GRB/hypernova remnants 
are promising candidates for the TeV unlDs. If TeV unlDs are 
hypernova remnants, we can predict more extended (and more 
numerous) TeV SNRs than observed, which may be discov- 
ered by expanding the search region to larger angular scales, 
even with the current instruments. We also propose a new 
process of TeV gamma-ray emission involving the decay of 
accelerated radioactive isotopes, such as 56 Co entrained by 
relativistic or semi-relativistic jets in GRBs/hypernovae. 
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